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TECHNICAL  NOTE  337^ 

TURBULENT -HEAT -TRANSFER  MEASUREMENTS 
AT  A  MACH  NUMBER  OF  2.0 6 
By  Maurice  J.  Brevoort  and.  Bernard  Rashis 

SUMMARY 


Turbulent-heat-transfer  measurements  were  obtained  through  the  use 
of  an  axially  symmetric  annular  nozzle  which  consists  of  an  inner  shaped 
center  body  and  an  outer  cylindrical  sleeve.  Measurements  taken  along 
the  outer  sleeve  gave  essentially  flat-plate  results  that  are  free  from 
wall  interference  and  comer  effects  for  a  Mach  number  of  2.06  and  for  a 

Reynolds  number  range  of  1.7  X  106  to  8.8  X  10?.  The  heat-transfer 
coefficients  are  in  good  agreement  with  available  data  from  V-2  r°cke  s 
and  also  check  the  Rubesin  theory  ana  the  Van  Driest  theory  for  a  Mac 
number  of  2.0  and  for  a  ratio  of  inner-surface  to  free -stream  temperature 
of  1  8.  The  temperature -recovery  factors  are  approximately  0.5  percen 
lower  than  the  factors  given  in  NACA  TN  2077  for  a  Mach  number  of  2.4. 

INTRODUCTION 


The  design  of  supersonic  aircraft  and  missiles  requires  engineering 
information  about  heat-transfer  coefficients  and  temperature -recovery 
factors  for  supersonic  speeds  that  extend  over  a  wide  range  of  Reynolds 
£  reference  lAocsl-heat-trsnsfer-coamcient 
presented  for  a  Mach  number  of  3.03-  Good  agreement  of  these  results  wibn 
theoretical  and  experimental  work  was  obtained.  This  method  of  testing 
and  reduction  of  Sta  is  readily  adaptable  to  obtaining  accurate  measure¬ 
ments  over  an  extended  range  of  both  Mach  and  Reynolds  numbers- 

The  purpose  of  this  investigation  is  to  extend  the  work  initiated  in 
reference  1  to  a  Mach  number  of  2.06.  The  same  type  of  apparatus  and 
method  S  reducing  the  data  vrere  used  in  this  investigation  as  vere  employed 
in  reference  1.  The  range  of  Reynolds  number  for  which  measurements  were 

obtained  is  from  1.7  X  10*  to  8.8  X  10T.  The  results  cover  a  tempera¬ 
ture  difference  of  approximately  20°  at  40  seconds  after  starting  to 
approximately  5°  at  120  seconds  after  starting.  The  average  value  of  the 
ratio  of  inner-surface  temperature  to  free-stream  temperature  Tw/Tco 
was  1.8. 
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SYMBOLS 


specific  heat  of  sleeve  material,  Btu/lb-°R 
specific  heat  of  air  at  constant  pressure,  Btu/lb-°R 


acceleration  due  to  gravity,  ft/sec^ 

heat -transfer  coefficient,  Btu/ft^-sec-°R 

heat  conductivity,  Btu/ft-sec-°R 

Mach  number 

Nusselt  number,  hx/k 

Prandtl  number,  pCpg/k 

Reynolds  number,  pVx/p 


Stanton  number,  . .  =  — £ — 

R  Pr  pVCpg 

average  wall  temperature,  °R 

effective  stream  air  temperature  at  wall,  some  tempera¬ 
ture  which  gives  a  thermal  potential  which  is  independent 
of  heat-transfer  coefficient  h,  °R 

stagnation  temperature,  °R 

inside-surface  temperature  of  nozzle  sleeve,  °R 
free -stream  temperature,  °R 
time,  sec 

free -stream  velocity,  ft/sec 

specific  weight  of  sleeve  material,  lb/sq.  ft 

longitudinal  distance  along  sleeve,  ft  (unless  indicated 
otherwise ) 
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p 

p 


T  -  T 

recovery  factor,  -§ - — 

T.  -  T 

t  °° 

dynamic  viscosity  coefficient,  lb-sec/sq  ft 
free-stream  density  of  air,  slugs /cu  ft 

APPARATUS  AND  METHOD 


The  apparatus  consisted  of  an  axially  symmetric  annular  nozzle  which 
was  directly  connected  to  the  settling  chamber  of  one  of  the  cold-air 
blowdown  jets  of  the  Langley  gas  dynamics  laboratory.  The  nozzle  had  a 
shaped  wooden  center  body  and  two  outer  sleeves.  The  first  sleeve  was 
constructed  from  8-inch-diameter,  extra  heavy,  seamless  carbon-steel  pipe 
and  the  second  was  constructed  from  l/l6-inch  stainless  steel  which  was 
rolled  into  a  cylinder  and  welded;  both  were  surfaced -machined  inside 
and  outside  to  a  wall  thickness  of  O.388  inch  and  0.060  inch,  respectively. 
A  detailed  drawing  of  the  apparatus  is  shown  in  figure  1  which  gives  the 
location  of  the  thermocouples  and  the  static -pressure  orifices.  Several 
static -pressure  orifices  and  thermocouples  were  located  90°  apart  around 
the  sleeve  to  calibrate  the  flow  for  axial  symmetry.  The  coordinates  of 
the  center  body  are  given  in  table  I.  The  center  body  was  designed  by 
using  three-dimensional  characteristics.  (See  ref.  2.)  Station  locations, 
shown  in  figure  1,  refer  to  the  distance  in  inches  from  the  nozzle  mini¬ 
mum  station. 

Details  of  the  installation  of  the  thermocouples  and  static-pressure 
orifices  are  shown  in  figures  2  and  3*  The  thermocouples  are  located 
0.060  inch  from  the  inner  surface  of  the  sleeve.  The  wires  are  JO-gage 
copper-constantan  (0.010  inch  in  diameter);  conduction  along  the  wire 
length  is  negligible.  As  indicated  in  figure  2,  the  thermocouples  are 
in  intimate  contact  with  the  steel  wall  so  that  thermal  resistance  at 
this  junction  is  negligible. 

The  Mach  number  distribution  was  measured  both  lengthwise  and  around 
the  sleeve.  The  orifices  located  around  the  sleeve  were  used  to  check  the 
alinement  of  the  center  body  within  the  sleeve.  The  results  are  shown  in 
figure  4  and  are  accurate  to  ±0.01. 

The  temperature -recording  equipment  consisted  of  three  synchronized 
high-speed  self-balancing  Brown  ElectroniK  strip  chart  recorders,  having 
a  total  of  38  switches,  two  of  which  were  connected  to  the  settling- 
chamber  thermocouples  and  the  others  were  connected  to  the  sleeve 
thermocouples . 
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A  typical  test  consisted  of  operating  the  nozzle  at  the  desired 
pressure.  The  stagnation  temperature  starts  at  essentially  room  tempera¬ 
ture  and  decreases  as  the  piping  is  cooled,  as  illustrated  in  figure  5 
where  stagnation  temperature  is  plotted  against  time  for  a  settling -chamber 
pressure  of  10 6  lb/sq  in.  gage.  The  wall  temperature  also  starts  at  essen¬ 
tially  room  temperature  and  tends  to  approach  the  equilibrium  temperature 
which  is  approximately  25°  F  below  stagnation  temperature.  This  variation 
is  shown  in  figure  6  where  wall  temperature  is  plotted  against  time  for 
stations  12  and  18  for  a  settling -chamber  pressure  of  106  lb/sq  in.  gage. 
Test  runs  were  made  for  settling -chamber  pressures  of  1,  106,  and 

185  lb/sq  in.  gage.  The  test  at  a  pressure  of  1  lb/sq  in.  gage  was  made 
with  a  sleeve  0.060  inch  thick,  and  the  test  setup  was  evacuated  to  an 
absolute  pressure  of  1.0  inch  of  mercury.  Except  for  the  first  20  seconds, 
the  pressures  were  maintained  constant  for  each  test  run.  The  recorders 
were  calibrated  immediately  before  and  after  each  run  and  were  found  to 
be  accurate  to  ±1°  F. 


In  figure  7  are  plotted,  for  various  times  during  the  test  run,  the 
values  of  wall  temperature  against  longitudinal  distance  along  the  cylinder. 
These  values  were  used  to  determine  the  rate  of  change  of  the  longitudinal 
d2 

conduction  k  — 2§X  along  the  cylinder.  Test  results  were  taken  only  for 
dx2 


d.2?T 

the  length  of  the  cylinder  for  which  k  - —  was  zero. 


dx£ 


REDUCTION  OF  DATA 


Short  test  runs  in  which  time  histories  of  the  wall  temperature  are 
obtained  give  the  essential  data  from  which  heat-transfer  coefficients 
and  temperature -recovery  factors  may  be  determined.  The  reduction  of 
these  data  requires  a  method  in  which  the  stagnation  temperature,  and 
accordingly  the  equilibrium  temperature,  may  vary  in  an  arbitrary  manner. 

The  Stanton  number  is  calculated  from 


St 


h 

pVcpg 


and  the  heat -transfer  coefficient  is  calculated  from 


(1) 


(2) 


NACA  TN  3374 


5 


By  definition,  the  recovery  factor  is 


^r 


-  T 

_ XL 


(3) 


The  method  of  reducing  the  data  is  simply  to  select  a  recovery  factor 
and  then  obtain  Tg  from  equation  (3).  Substitute  this  value  of  Te  into 

equation  (2)  and  obtain  values  of  h  for  different  heat-flow  rates.  These 
values  of  h  are  then  substituted  into  equation  (l).  The  true  values  of 
Tg,  T]r,  and  St  are  obtained  when  St  is  constant  with  time  (for  dif¬ 
ferent  heat-flow  rates ) .  Figure  8  shows  the  values  used  in  evaluating  the 
Stanton  number  and  recovery  factor  at  station  12  for  a  settling-chamber 
pressure  of  106  lb/sq  in.  gage. 

The  values  of  specific  heat  and  specific  weight  of  the  sleeve  mate¬ 
rial  were  taken  from  reference  3.  The  values  of  viscosity  and  Prandtl  num¬ 
ber  (O.7I)  for  air  were  also  taken  from  reference  3  and  were  based  upon 
the  inside -surf ace  temperature  of  the  sleeve.  The  value  of  Tw  chosen 

was  at  80  seconds  after  starting.  (This  arbitrary  choice  of  temperature 
can  be  made  since  the  viscosity  coefficient  is  practically  a  constant 
factor  of  the  heat -conduction  coefficient.) 


PRECISION  OF  EQUIPMENT 


The  basic  data  consist  of  time  histories  of  temperature  at  each 
thermocouple  location.  The  recorders  were  calibrated  immediately  before 
and  after  each  run  and  were  found  to  be  accurate  to  ±1°  F.  The  thermo¬ 
couple  wires  are  30-gage  (0. 010 -inch-diameter)  copper -constantan  and  of 
such  small  diameter  that  conduction  along  their  length  can  be  neglected. 
The  wires  axe  insulated  except  at  the  junction  where  intimate  contact 
with  the  sleeve  keeps  any  thermal  resistance  to  a  minimum. 

The  sleeve  is  enclosed  in  a  vacuum  jacket  to  avoid  free  convection 
and,  for  the  temperature  range  considered,  radiation  effects  range  up  to 
2  percent.  The  product  of  the  specific  weight  and  the  specific  heat  of 
the  sleeve  material  as  used  in  equation  (2)  is  accurate  to  ±5  percent. 

cLT 

The  value  of  used  in  equation  (2)  is  that  of  the  actual 

dt 

thermocouple  reading,  whereas  it  should  be  the  value  associated  with  the 
average  temperature  in  the  thickness  of  the  sleeve  at  each  point.  During 
the  test  runs  the  inner-surface  temperatures  are  lower  than  the  average 
temperatures  by  approximately  1°  to  I.50  at  the  beginning  of  the  runs 
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(40  seconds  after  starting)  and  become  the  same  at  the  end  of  the  runs 
(150  seconds).  This  deviation  causes  the  derivative  of  inner-surface 
temperature  with  time  to  be  less  than  the  derivative  of  average  tempera¬ 
ture  with  time  by  approximately  1  percent. 

The  analysis  of  reference  4  was  used  to  determine  the  effect  of  the 
sleeve  diameter  upon  the  flow  characteristics .  This  effect  was  computed 
to  be  less  than  one-hundredth  of  1  percent. 


RESULTS  AND  DISCUSSION 


Figure  7  shows  the  variation  of  wall  temperature  with  longitudinal 
distance  along  the  cylinder  for  a  settling -chamber  pressure  of 
106  lb/sq  in.  gage.  Over  the  test  range,  the  wall  temperatures  are  con¬ 
stant  along  x.  In  the  relation 

vc  dTav/dt 
Tw  -  Te 

w  and  c  are  constant  and  dTav/ dt  is  constant  because  Tw  is  constant. 

Therefore,  if  there  is  to  be  a  variation  of  h  with  Reynolds  number  or 
x,  Te  must  vary  with  x.  The  value  of  Te  obtained  for  the  test  at  a 

settling-chamber  pressure  of  106  lb/sq  in.  gage,  evaluated  at  80  seconds 
after  starting,  actually  decreases  approximately  1.5°» 

Figure  9  shows  the  variation  of  local  Nusselt  number  with  local 
Reynolds  number.  The  value  of  x  used  in  evaluating  these  numbers  was 
based  upon  an  x  that  was  considered  to  be  zero  at  the  nozzle  minimum 
station.  A  single  line  faired  through  all  the  data  points  has  a  slope  of 
approximately  0.8;  but  if  lines  are  faired  through  the  points  for  each 
settling-chamber  pressure,  the  slopes  are  slightly  higher.  It  is  reason¬ 
able  to  assume,  however,  that  the  higher  stagnation  pressures  would  pro¬ 
duce  earlier  transition  and  the  zero  value  of  x  would  move  upstream 
with  increasing  stagnation  pressure.  Hence,  if  it  is  assumed  that  data 
from  each  settling-chamber  pressure  should  merge  into  a  continuous  line 
(rather  than  a  series  of  steps),  x  may  be  adjusted  for  zero  x-locations 
(the  effective  beginning  of  the  turbulent  boundary  layer)  which  would 
make  all  the  test  runs  coincident.  This  adjustment  has  been  made  in 
figure  10  by  using  values  of  x  equal  to  zero  at  1.60  inches  downstream 
of  the  nozzle  minimum  station  for  a  settling-chamber  pressure  of 
1  lb/sq  in.  gage  and  2.00,  4.40,  and  6.00  inches  upstream  of  the  minimum 
for  the  chamber  pressures  of  38,  106,  and  185  lb/sq  in.  gage,  respectively. 
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The  Nusselt  numbers  were  found  to  vary  from  1,690  to  44,950  for  the 

Reynolds  number  range  of  1.7  x  10^  to  8.8  x  10^  (based  on  adjusted 
zero  x-locations ) .  For  comparison,  the  curves  based  upon  the  analyses 
of  Van  Driest  (ref.  5)  and  Rubesin  (ref.  6)  are  shown.  In  the  Van  Driest 
analysis  was  considered  to  be  1.8  (the  average  value  of  the  test 

results).  For  comparison  with  these  references,  the  data  were  computed 
by  using  free -stream  temperature  to  determine  the  density  and  the  velocity. 
The  wall  temperature  was  used  to  determine  the  viscosity  and  the  Prandtl 
number.  The  experimental  values  from  the  data  of  V-2  rockets  (ref.  7) 
are  also  determined  in  this  manner.  The  Mach  number  range  for  the  V-2  data 
in  figure  10  is  from  1.95  1°  2.07«  The  agreement  between  the  results  of 
this  investigation  and  those  of  references  5,  6,  and  7  is  good. 

Figure  11  shows  the  variation  of  local  temperature -recovery  factor 

with  local  Reynolds  number.  The  variation  is  from  O.89I  at  R  =  1.7  x  10^ 

to  O.87O  at  R  =  8.8  x  10?.  The  results  are  compared  with  a  curve  which 
represents  the  fairing  of  the  data  of  reference  8.  The  data  of  reference  8 
are  for  a  Mach  number  of  2.4.  The  results  obtained  in  this  investigation 
are  approximately  0.5  percent  lower.  Also  included  for  comparison  are  the 

curves  for  recovery  factor  equal  to  Pr1/?  and  Pr1/2.  The  wall  tempera¬ 
ture  was  used  to  determine  the  Prandtl  number. 


CONCLUDING  REMARKS 


Turbulent-heat-transfer  measurements  that  gave  essentially  flat -plate 
results  were  obtained  for  a  Mach  number  of  2.06  and  for  a  Reynolds  number 

range  of  1.7  x  10^  to  8.8  x  10^.  The  Nusselt  numbers  are  in  good 
agreement  with  theoretical  analyses  and  with  experimental  data  obtained 
with  the  V-2  rockets.  The  temperature -recovery  factors  are  approximately 
0.5  percent  lower  than  other  available  data  for  a  Mach  number  of  2.4. 


Langley  Aeronautical  laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va. ,  December  3,  1954. 
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Four  thermocouples  at  stations 


1 


MCA  TN  3374 


Figure  10.-  Variation  of  local  Nusselt  number  with  local  Reynolds 
number  for  corrected  locations  of  x  =  0.  Viscosity  and  Prandtl 
number  determined  for  wall  temperatures . 
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